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ABSTRACT 

Various  efficient  algorithms  for  detecting  intersections  and  computing 
closest  neighbors  in  a  set  of  convex  bodies,  are  presented  and  analyzed.  These 
algorithms  generalize  known  techniques  for  solving  these  problems  for  sets  of 
points  or  of  line  segments.  Of  particular  interest  is  the  construction  of  a  general- 
ized Voronoi  diagram  for  a  set  of  n  circular  bodies  in  time  0{n  log"n),  and  its 
applications. 


1.   Introduction 

Let  S  be  a  set  of  n  closed  convex  two-dimensional  bodies  of  relatively  simple  structure  (eg 
circular  discs,  straight  segments,  polygons  of  few  sides,  or  expansions  by  some  amount  of  such 
objects).  In  thb  paper  we  consider  a  variety  of  problems  associated  with  such  a  set  5.  Typical 
such  problems  are: 

I.  Do  any  two  objects  in  S  intersect? 

II.  More  generally,  suppose  that  we  assign  a  'color'  to  each  object  in  S.  Does  there  exist  a  pair  of 
objects  in  5  having  different  colors  and  intersecting  each  other? 

III.  If  no  two  objects  in  5  intersect,  what  is  the  smallest  distance  between  any  two  objects  in  5 
(or,  in  the  'colored'  version,  what  is  the  smallest  distance  between  two  objects  in  S  having 
different  colors)?  More  generally,  for  each  object  B  in  S  find  the  object  in  S  (of  different  color) 
nearest  to  B. 

rV.  Preprocess  5  so  that,  given  an  arbitrary  'query  point"  X,  the  object  in  5  nearest  to  X  can  be 
found  quickly. 

Problems  of  this  sort  arise  in  robotics  and  are  related  to  the  problem  of  detecting  and  avoid- 
ing collisions  between  a  moving  subpart  of  a  robot  system  and  stationary  objects,  or  between  two 
or  more  moving  subparts  of  such  a  system.  In  this  note  we  will  simplify  the  problem  by  assuming 
that  each  of  the  robot  subparts  and  the  stationary  obstacles  is  either  a  closed  convex  object  of  a 
simple  form,  or  else  is  covered  by  finitely  many  such  objects.  Note  that  in  the  'colored'  setting 
some  of  the  objects  involved  in  our  problem  may  be  known  a  priori  to  intersect  (e.g.  circles  cover- 
ing two  robot  subparts  hinged  together,  or  two  objects  covering  some  robot  subpart  which  overlap 

'*'  This  research  has  been  supported  in  part  by  CNR  Grant  N00014-82-K-0381  and  by  a  grant  from  the  U.S.- 
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each  other  may  always  intersect).  The  'colored'  version  of  our  problems  allows  for  this  situation 
by  looking  only  for  intersection  of  subparts  which  would  not  intersect  under  normal  conditions 
(e.g.  subparts  belonging  to  two  distinct  robot  'arms',  or  a  robot  subpart  and  an  obstacle,  etc.). 

Efficient  solution  of  these  problems  in  the  three-dimensional  case  would  facilitate  construc- 
tion of  an  'off-line'  debugging  system  for  robot  control  programs  to  check  whether  collision  occurs 
along  a  planned  path  of  motion,  and  would  also  make  it  more  feasible  to  check  in  real  time 
whether  a  moving  subpart  of  the  system  is  getting  dacioicusly  close  to  another  (moving  or  sta- 
tionary) object.  The  study  of  the  two-dimensional  case  cirritd  out  in  this  paper  is  intended  to 
suggest  ideas  which  may  be  useful  in  attacking  tL«r  three-dimensional  case. 

2.   Detecting  Intersection  of  Convex  Obji;jis 

In  this  section  we  consider  the  first  problem  posed  above  and  show  that  it  can  be  solved  in 
time   0(n  log  n)  by  a  straightforward  modification  of  '<iZ.  algorithm  due  to  Shamos  [Sh]  which 
tests  for  intersection  of  straight  line  segments.    We  assume  thac  I  le  bi>dies  in  5  are  convex  and 
have  simple  enough  structure  so  that  each  of  the  foUowiag  opetatsoat,  takes  constant  time: 
(i)  Check  whether  two  specific  bodies  in  S  intersect  each  other, 
(ii)  For  each  B  e  S,  find  the  smnllest  ard  largest  abscisaaiB  of  point*  in  B. 
(iii)  For  each  B  £  S,  find,  for  a  given  .ibscma  x,  a  poiiit  (l,j/)  G  B. 

Let  Bi,  .  .  .  ,B„  be  the  objects  i;^  S.  For  eachV=i,  ■  •  •  . "  'e'  "j  (r«9P-  *;)  ^e  the  smallest 
(resp.  largest)  abscissa  of  a  point  in  B,.  For  suniplicity  we  assaine  that  the  2n  numbers 
aj,  bj,  ;  =  1,  .  .  .  ,  n,  are  all  distinct  (see  Fig.  1).-  ^-  '"       '''P  '•-•'• 


L(x,)  =[83.82.8,] 


Fig.  1.  An  instance  of  the  intersection  problem. 

The  main  idea  used  by  the  algorithm  described  below  is  that  if  we  draw  a  vertical  line 
z  =  zc.  it  will  cut  (some  of)  the  objects  in  S  in  straight  segments  which,  if  the  objects  do  not  inter- 
sect each  other  along  that  line,  are  disjoint  from  each  other.  Hence  the  objects  in  5  which  inter- 
sect the  line  z=Xo  c=^  ^e  linearly  ordered  in  a  list  L(io)  '^  which  an  object  5,  precedes  another 
B,  if  the  segment  cut  off  S,  by  x  =  Xq  lies  below  the  corresponding  segment  cut  off  Bj.  Note  that 
since  the  objects  in  5  are  convex,  the  list  L(xq)  remains  unchanged  as  z^  increases  until  either  the 


line  x=Xo  meets  a  new  object  Bt  (this  will  happen  when  io=at),  or  when  it  stops  making  contact 
with  some  object  B^  (just  after  Xo=6t),  or  when  two  of  the  objects  in  L{xo)  intersect  at  x=xq. 
Moreover,  the  leftmost  intersection  (if  any  exists)  of  any  two  objects  B,  B'  6  5  will  occur  at 
some  z  =  io  such  that,  for  i  slightly  less  than  Iq,  the  list  L(X)  contains  B  and  B'  as  adjacent  ele- 
ments. (Here  we  ignore  the  special  case  in  which  B  and  B'  intersect  at  a  point  which  is  leftmost 
in  one  of  these  objects,  which  requires  a  slightly  different  argument.) 

In  view  of  these  observations,  to  detect  the  presence  or  absence  of  an  intersection  one  sim- 
ply has  to  check  repeatedly  whether  any  two  adjacent  elements  in  the  list  L(xq)  intersect  each 
other.  Plainly,  these  checks  have  to  be  performed  only  at  points  where  L{xq)  changes  (i.e.  at  the 
points  XQ=aj,  bj,  j  =  l,  .  .  .  ,  n),  and  one  only  needs  to  test  newly  adjacent  pairs  in  L(xq).  To 
facilitate  execution  of  these  operations,  the  list  L  ?an  be  maintained  as  a  2-3  tree,  allowing  all  the 
required  list-main  tain  ance  operations  to  be  performed  in  time  0(log  n)  (see  below  for  details). 

Details  are  as  follows.  The  algorithm  begins  by  sorting  the  2n  numbers  aj,bjj  =  l,  .  .  .  ,n, 
in  increasing  order,  and  then  processes- them  from  left  to  right.  Initially,  the  list  L  is  empty. 
Suppose  that  the  abscissa  curreptly  being  processed  is  one  of  the  a^.  Then  L  is  updated  by 
inserting  the  object  Bj  into  L  in  its  proper  place,  using  a  standard  2-3  tree  search  during  which 
comparison  of  two  objects  B,  B'  is  acco'mpnshed  by  comparing  two  representative  points  in  the 
intersection  of  x  =  aj  with  B;  Bf  respectively  (we  have  assumed  that  this  can  be  done  in  constant 
time).  After  insertion,  the  algorithm- |nds  the  two,  objects  B,  B  immediately  preceding  and 
succeeding  Bj  \a  L  ,  and  checks  whether  either  B  or  B'-  intersects  Bj. 

Similarly,  if  the  abscissa  curr.ently,  being  processed  is  6^,  then  the  object  Bj  is  deleted  from 
L,  using  essentially  the  technique  just  outlined.  After  deletion,  the  algorithm  finds  the  two 
objects  in  L  which  immediately  preceded  and  followed  B,  prior  to  its  deletion  (these  will  have 
become  newly  adjacent  in  L  after  deletion  of  B,),  and  determines  whether  they  intersect  each 
other. 

The  algorithm  halts  whenever  an  intersection  is  detected,  or,  if  no  intersection  has  been 
detected,  when  all  the  abscissae  a,  and  bj  have  been  processed,  in  which  case  the  algorithm 
reports  that  there  is  no  intersection  between  the  objects  in  S. 

The  correctness  of  the  algorithm  follows  from  the  preceding  observations.  The  time  com- 
plexity of  the  algorithm  is  0(n  log  n)  since  processing  of  each  of  the  2n  abscissae  a.,bj  can  be 
accomplished  in  0(log  n)  time,  using  a  2-3  tree  representation  for  the  list  L  . 

3.   Voronoi  Diagrams  for  Circular  Objects 

The  algorithm  presented  in  the  preceding  section  does  not  solve  the  more  complicated 
'colored  object'  intersection  problem  posed  in  the  introduction.  Indeed,  the  argument  justifying 
the  correctness  of  the  algorithm  breaks  down  as  soon  as  an  intersection  is  detected,  so  that  if  the 
first  intersection  detected  is  between  two  objects  having  the  same  color  we  can  no  longer  use  the 
procedure  described  to  find  additional  intersection  points.  To  handle  the  colored  intersections 
problem,  we  therefore  propose  a  different  approach  based  on  generalized  Voronoi  diagrams.  We 
will  show  such  an  approach  which  can  be  used  to  handle  the  special  case  where  all  the  objects  in 
the  set  5  are  circular  discs,  not  necessarily  of  equal  radii. 

Let  each  of  the  objects  S^  6  5  be  a  disc  of  radius  r,  about  the  center  z,,  for  j  =  l,  .  .  .  ,n. 
These  circular  discs  need  not  be  disjoint  from  each  other,  and  may  intersect,  or  even  contain,  one 
another.  We  define  a  generalized  Voronoi  diagram  Voro(5)  associated  with  the  set  5  as  follows. 
For  each  iy^j  define 

H{iJ)={y  e  ET-.  d(x„y}-r,<d(x„yy-r^}, 

i.e.  the  set  of  all  points  whose  distance  from  B,  is  no  greater  than  their  distance  from  B^  (note 
that  the  dbtance  of  y  from  B,   is  taken  as  the  distance  of  y  from  the  boundary  of  B,  with  a 


positive  sign  if  y.  lies  outside  B,  and  with  a  negative  sign  if  y  lies  inside  B,).  Then  define  the 
(closed)  Voronoi  cell  V{i)  associated  v/ith  B,  to  be 

i.e.  the  set  of  all  points  y  whose  distance  from  B,  'a  no  ^-eater  than  y'a  distance  from  any  other 
element  of  S.  (We  will  sometimes  refer  to  the  point  2,  as  the  center  of  this  Voronoi  cell;  hence 
the  center  of  V'(»)  is  the  same  as  the  center  cf  the  dbc  B,.)  Finally,  the  Voronoi  diagram  VorQ{S) 
is  defined  to  be  the  set  of  points  which  belong  to  more  than  one  Voronoi  cell.  For  simplicity  we 
assume  that  no  point  in  Voro(5)  lies  in  more  than  four  Voronoi  cells.  This  assumption  generalizes 
the  familiar  assumption  concerning  Voronoi  dingrams  associited  with  a  set  of  points,  in  which  one 
requires  that  no  more  than  four  of  these  points  be  cocircular.  Fig.  2.  shows  an  e.xample  of  such  a 
Voronoi  diagram. 


V(l) 


Fig.  2.  The  Voronoi  diagram  of  a  set  of  circular  discs. 

Generalized  Voronoi  diagrams  have  been  previously  introduced  arA  analyzed  by  various 
authors;  see  e.g.  Kirkpatrick  [Ki|,  who  considers  Voronoi  diagrams  for  a  set  of  line  segments  and 
points,  and  Lee  and  Drysdale  [LD|,  who  consider  Voronoi  diagrams  for  sets  of  line  segments  and 
(or  set"  of  circular  objects.  Our  setup  is  somewhat  different  from  that  of  [LDj,  in  that  we  allow 
the  objects  in  S  to  intersect  each  other,  or  even  to  contain  one  another,  whereas  in  [LD]  the  cir- 
cles are  required  to  be  disjoint.  Lee  and  Drysdale  present  an  0(n  log'n)  algorithm  for  construct- 
ing the  Voronoi  diagram  of  n  disjoint  circles,  but  their  algorithm  is  given  in  a  very  sketchy  form 
without  any  analysis  of  the  structure  of  the  diagram,  and  without  any  proof  of  its  correctness.  In 
addition,  they  assume  availability  of  an  0(r.  log  n)  algorithm  for  locating  n  points  in  general 
(nonlinear)  planar  subdivisions,  but  give  no  details  concerning  such  an  algorithm  (the  algorithm  of 
Preparata  [Pr]  lo  which  they  refer  has  to  be  modified  to  handle  nonlinear  arcs,  and  is  unneces- 
sarily complicated  for  the  task  of  simultaneous  location  of  n  points).  A  by-product  of  the  algo- 
rithm that  we  present  in  Section  4,  is  such  a  relatively  simple  algorithm,  which  may  also  be  of 
independent  interest. 

The  generalized  Voronoi  diagram  just  defined  has  the  following  properties. 
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(1)  The  collection  of  Voronoi  cells  covers  the  whole  plane. 

Proof:  Immediate.  Indeed,  given  any  y  G  £^  ,  it  will  belong  to  the  ceil  V{i)  for  which 
<'(^i.v)-'"i=in'n{('(3;;,I/)-r;  :  j  =  l,  .  .  .  ,n}. 

(2)  V(i)  is  empty  iff  B,  is  wholly  contained  in  the  interior  of  another  disc  B/,  V(i)  has  an  empty 
interior  iff  B,  is  wholly  contained  in  another  disc  By 

Proof:  For  the  first  assertion,  note  that  if  V(i)  is  not  empty  it  contains  a  point  x  such  that 

d(x,x,}-r,  <  min{d{^x,Xj)-rj  :  j  =  l,  .  .  .  ,«}. 

The  triangle  inequality  then  show»tiiat  this  same  ioequality  holds  for  x  =  x,,  i.e. 

-r,=m\n{d{x,,Xj)-rj  :  j==l,  .  .  .  ,n}, 

which  implies  that,  for  any  of  the  other  circles  Bj  of  5,  the  point  of  B,  at  ma.ximum  distance 
from  Xj  is  not  interior  to  B,.  Thus  there  is  no  B,  whose  interior  contains  fl, .  Conversely,  if  V(i) 
is  empty,  then  x,  ^  V(i),  so  that  there  exists  j^i  such  that 

i.e. 

d(x,,x,)+r,<rj 

which  is  to  say  that  Bj  contains  B,  in  its  interior.  Similarly,  if  the  interior  of  V(i)  is  not  empty, 
it  contains  a  point  z  such  that 

d(x,x,)-r,  <  mia{d(x,Xj)-rj  :  ;'=1,  .  .  .  ,n}  , 

and  again  this  inequality  must  hold  for  x=x,.  Thus 

-r,<d(x,,Xj)-rj,         for  each  ;■  7^1, 

Hence  rj<d(x,,Xj)+  r,  for  each  ;,  i.e.  the  point  of  B,  at  maximum  distance  from  x,  does  not 
belong  to  B,.   Thus  there  is  no  B,  which  contains  B,. 

Conversely,  if  the  interior  of  V(i)  is  empty,  x,  does  not  belong  to  this  interior,  so  that  there 
exists  jy^i  such  that 

d(x„x,hr,<-r, 
or 

d(x„x,)<r,-r„ 
which  is  to  say,  B,  is  wholly  contained  in  B,. 

(3)  Vor^S)  consists  of  straight  or  hyperbolic  arcs. 

Proof:  Let  y  6  V'oro(5)  be  such  that  y  lies  in  both  Voronoi  cells  V(i)  and  V(j).  Then  we  have 

d{x,,y)-d(x,,y)=r,-r, 

The  locus  of  points  satisfying  this  condition  is  a  hyperbolic  arc  having  x,  and  i.  as  foci,  or,  if 
r,=rj,  the  perpendicular  bisector  to  the  segment  [i,,!;].  (Note  also  that  this  (generally  hyper- 
bolic) locus  degenerates  into  a  half-line  if  d(x,,Xj)=±  (r.-r^),  i.e.  if  one  of  the  discs  S,,  B,  wholly 
contains  the  other. 

(4)  Each  nonempty  Voronoi  cell  V(i)  is  star-shaped  with  respect  to  the  point  j, .  Moreover,  if  V{i) 
has  nonempty  interior,  then  the  interior  of  a  segment  connecting  x,  to  a  point  on  the  boundary  of 
V(i)  does  not  intersect  the  interior  of  any  other  Voronoi  cell,  and  such  a  segment  can  intersect 
another  Voronoi  cell  V(j)  only  if  the  corresponding  disc  Bj  is  wholly  contained  in  B,  (so  that,  by 


remark  (3),  V(j)  has  empty  interior). 

Proof:  Let  y  6  V'(j),  and  let  /  be  the  segment  connecting  j,  to  y.  We  first  claim  that  each  point 
z  in  the  interior  of  /  is  contained  in  V(i).  Indeed,  if  this  were  false  then  there  would  exist  a  point 
;  in  the  interior  of  /  which  does  not  belong  to  V'{»).  By  (1)  there  would  exist  j^i  such  that  z 
belongs  to  V(j).  Hence 

d(xj,z)-r.<d(x,,z}-r,. 
By  the  triangle  inequality  (and  since  ;  lies  betwcca  s,  zxid  y  on  &  straight  line)  we  have 

d(x,,y}-r,<d{x^,z)+  d(z,y)-r.<d{x,,z]-]-  diz ,y)-r.  =  d{x,,y}-r,. 
Thus  y  cannot  belong  to  V(i),  a  contradictioL  Wiiich  proves;  th;it  V{i)  is  star-shaped. 

Next  suppose  ihat  V{i)  has  nonempty  interior.  Let  y  €  V'yt)  and  I  be  as  above,  and  sup- 
pose that  /  contains  an  interior  point  z  wnich  also  belongs  to  some  other  Voronoi  cell  V(j).  Then 
we  have 

d{xj,z)-i'j  =  d{x,,2)-r,. 
Using  the  triangle  inequality  as  before,  we  obtain 

d{xj,y)-r.<d(.v^,r)+  (i(z,i/)-r^=f=^^g,,,^);t-  d(z.y)-r, 

-^d[x,,y}-r,<dix,)i)-rj, 

since  y  G  V(i).  Hence  ij  lies  on  the  line  contrining  I, ,aR^,^^Iies  between  i,  and  y  on  this  line. 
However,  this  implies  that  ,.      ',!,/,, 

d(x„x,)==±(r,-r^l    ■      ■- 

i.e.  that  one  of  the  discs  B,,B,  contains  wholly  the  other.  But  by  (2)  B  is  not  wholly  contained 
in  any  other  disc,  so  that  Bj  is  wholly  contained  in  B,,  and  consequently  V{j)  has  empty  interior. 
This  establishes  our  two  final  assertions. 

(5)  The  intersection  /  of  three  Voronoi  cells  V'(»),  V(j),  V(k)  can  consist  of  at  most  finitely  many 
points;  moreover,  if  all  these  three  cells  have  nonempty  interiors,  then  /  consists  of  at  most  two 
points. 

Proof:  We  prove  only  the  second  assertion,  which  is  needed  below,  since  the  first  assertion  is 
much  simpler.  Let  y  be  a  point  in  /.  By  (4)  the  interior  of  the  segment  L,  (resp.  Lj.L^  )  con- 
necting y  with  J,  (resp.  x,,Xi;)  is  contained  in  V(»)  (resp.  V(j),V(k))  and  in  no  other  cell.  Let  y' 
be  another  point  in  /,  and  let  L,'  ,Lj'  ,Li['  be  the  segments  connecting  y'  with  j,,  Xj,  j^  respec- 
tively. It  is  clear  that  no  two  of  these  six  segments  can  intersect  one  another  (except  at  an  end- 
point).  It  follows  that  /  cannot  contain  a  third  point  z,  because  at  least  one  of  the  segments  con- 
necting z  to  the  three  centers  x^,Xj,Xi  would  have  to  intersect  one  of  the  preceding  six  segments, 
which  is  impossible. 

Definition:  The  modified  Voronoi  diagram  Vor(S)  is  defined  to  consist  of  the  boundaries  of  all 
cells  in  Votq(S)  having  nonompty  interior.  All  other  ceils  are  discarded  from  the  modified 
diagram. 

(6)  Let  D'  denote  the  dual  of  Vor[S),  defined  to  be  a  graph  whose  vertices  are  the  points 
X,,  1  =  1,  .  .  .  ,n  for  which  V'(«)  has  nonempty  interior,  and  which  contains  an  edge  [x,,2j  if  V'(j) 
and  V[j)  have  a  nonempty  intersection  including  at  least  one  open  arc.  If  V(i)  and  V[j)  intersect 
in  more  than  one  one  arc,  define  the  graph  D'  to  contain  multiple  edges  connecting  z,  to  x,,  one 
for  each  such  arc.  Then  D'  is  a  planar  graph,  each  of  whose  faces  (except  for  the  outer  one) 
contains  at  least  three  edges. 

Proof:  We  define  an  embedding  ot  D'  in  £^  as  follows.  Each  vertex  7,  is  mapped  to  itself  (as  a 
point  in  E^  ).  Let  e  =  [2,,rj  6  £*'    be  an  edge  corresponding  to  an  open  arc  ci  in  the  intersection 


of  V(i)  and  V(j).  It  follows  from  (5)  that  there  must  exist  at  least  one  point  y  £  a  which  does 
not  belong  to  any  other  Voronoi  cell.  We  then  map  the  edge  e  to  the  path  consisting  of  the  two 
segments  \x,,y]  and  \y ,x,].  To  see  that  the  resulting  graph  G  is  indeed  a  planar  embedding  of  D', 
suppose  to  the  contrary  that  two  distinct  edges  ei=\x,,Xj\,  e2=[it.i(|  ot  D'  map  to  paths  which 
intersect  each  other  at  a  point  z  which  is  not  a  vertex  of  G.  Let  z^  (resp.  z^)  be  a  point  on  the 
common  boundary  of  V(i)  and  V(j)  (resp.  V{k)  and  V(l))  such  that  e,  (resp.  e^  )  appears  in  G 
as  the  union  of  the  segments  [x,,Zi\  and  [zi,x,\  (resp.  \xy,,z^  and  (^i.^/!)-  Suppose  without  loss  of 
generality  that  the  segments  »=[z,,Zi]  and,  s\=^[xt,z^  intersect  at  a  point  other  than  a  common 
endpoint  i,  =Jt.  Since  by  (4)  the  interior  of  the  segment  s  (resp.  s')  is  contained  in  V{i)  (resp. 
V{k))  and  in  no  other  Voronoi  cell  having  nonempty  interior,  we  must  have  either  i  =  k  or  ^1=^2 
.  In  the  first  case  8  and  9'  meet  at  z,,  and  so  if  they  meet  at  another  point  they  must  overlap 
each  other,  which  is  possible  only  if  either  Zi=Z2  or  if  one  of  these  points  (say  Zi)  lies  in  Che  inte- 
rior of  the  other  segment  «'.  The  latter  assumption  would  contradict  the  fact  that  the  interior  of 
3'  b  wholly  contained  in  V(i)  and  in  no  other  cell  having  nonempty  interior.  Thus  we  must  have 
Zi  =  jo,  and  then  by  the  choice  of  these  points  it  follows  that  ;==/  too.  Then  it  b  plain  that  the 
two  arcs  of  the  common  boundary  of  "y(i)  and-  V^)  which  define  our  two  paths  are  identical. 
Hence  the  two  edges  ej  and  cj  are  not  distinct,  contrary  to  assumption.  All  this  shows  that  G  is  a 
planar  embedding  of  ZJ  *. 

Next  let  ^  be  a  face  of  G  "boondfefl^y  just  t\^o  "edges.  Then  these  edges  must  connect  the 
same  two  points  x,  and  x^,  through  two>  poipts -21,2^  eac|i  lying  on  a  different  open  arc  Oj,  ao  of 
the  common  boundary  of  V(i)  and  V(j).  But  then  Q  is  a  quadrangle  whose  vertices  are 
x,,Zi,Xj  and  Zn,  and  must  contain  the  point  v  at  which  one  of  the  arcs  aj,  Oo  ends,  which  must  be 
a  Voronoi  vertex.  Clearly  v  must  belong  to  the  three  celb  V{i),  V(j)  and  V(k),  where  5^  b  some 
other  disc  than  B,  and  Bj.  It  follows  that  x^  abo  lies  inside  Q  (because  the  interior  of  the  seg- 
ment connecting  x^  to  t;  lies  wholly  in  the  interior  of  V{k),  and  because  the  sides  of  Q  are  wholly 
contained  in  the  union  of  V(i)  and  V(y)).  Moreover,  since  we  have  assumed  that  no  Voronoi  ver- 
tex belongs  to  more  than  three  cells,  it  follows  that  y  is  a  vertex  at  which  three  Voronoi  edges 
meet,  one  of  which  separates  V(i)  from  V{k)  and  another  separates  V{k)  from  V{j).  Hence  G 
contains  two  paths,  which  are  images  of  edges  \x,,Xt\  and  [it.^J.  and  which  are  abo  contained  in 
Q.  This  implies  that  Q  b  not  a  face  of  G,  contrary  to  assumption,  establishing  the  second  asser- 
tion made  above. 

(7)  Since  D'  contains  0(n)  vertices,  and  since  each  of  its  faces  contains  at  least  three  edges,  it 
follows  by  Euler's  formula  that  it  has  at  most  0{n)  edges,  that  is,  Vor(S)  consists  of  at  most 
0{n)  connected  straight  or  hyperbolic  arcs. 

(8)  Let  C  be  the  convex  hull  of  the  union  of  all  the  discs  B,  G  S.  Note  that  the  boundary  of  C 
consists  of  an  alternating  sequence  of  straight  segments  and  circular  arcs,  the  circular  arcs  being 
boundary  portions  of  some  of  the  discs,  whereas  the  straight  segments  are  tangents  to  a  pair  of 
dbcs  in  5.  We  will  say  that  S,  and  B,  are  adjacent  along  the  boundary  of  C,  if  thb  boundary 
contains  a  straight  segment  tangent  to  both  B,  and  Bj.  Then  the  unbounded  edges  of  Vor{S)  are 
those  edges  that  are  common  to  two  celb  V(»),  V(j)  for  which  B,  and  Bj  are  adjacent  along  the 
boundary  of  C. 

Proof:  Let  e  be  an  unbounded  edge  of  C,  common  to  two  Voronoi  cells  V{i)  and  V(j).  Since  e 
is  either  a  straight  or  hyperbolic  arc,  it  tends  asymptotically  to  some  half-line  /.  Suppose,  without 
loss  of  generality,  that  /  is  the  positive  j/-axis,  and  let  a  =[0,^1  be  a  point  on  /.  For  sufficiently 
large  t  and  for  any  k  =  l,  .  .  .  ,n,  d{a,Bi]  behaves  asymptotically  as  i-»7t-ri,  where  r/j  is  the  y- 
coordinate  of  Xi,,  i.e.,  as  d(a,Btt)  where  5^/  b  the  image  of  B,  translated  parallel  to  the  i-axis 
until  its  center  lies  on  the  jz-axis.  It  follows  that 

for  every  k=l,  .  .  .  ,n.  Thb  however  is  easily  seen  to  imply  that  B,  and  Bj  are  adjacent  to  each 


other  along  C,  since  the  line  i/=»Ji+  r.  is  tangent  to  both  discs  B,,B,  and  all  the  discs  in  5  lie  in 
the  lower  half-plane  which  this  line  bounds,  so  that  the  portion  of  this  line  between  its  points  of 
tangency  with  S,  and  B,  belongs  to  the  boundary  of  C .  (An  extreme  case  that  we  need  consider 
is  that  in  which  the  line  y  =  ri,+  r,  is  also  tangent  to  a  third  disc  Bt,  at  a  point  lying  between  its 
points  of  contact  with  B,  and  B,.  However,  if  such  a  situation  arises,  it  is  easy  to  check  that  all 
the  points  on  e  sufficiently  far  away  are  nearer  to  St  than  to  one  of  S,,  B,,  contradicting  the 
definition  of  e.  The  converse  statement,  namely  thjt  any  paiv  cf  adjacent  discs  along  C  induces 
an  unbounded  Voronoi  edge,  can  be  proved  using  the  above  argument  in  reverse. 

(9)  Vor(S)  need  not  be  connected.  In  fact,  it  r^a'i  'lav*;  un  i;o  0{n)  connected  components.  How- 
ever, every  connected  component  of  Vot(S)  is  unbounded. 

Proof:  Cocsider  the  following  set  5  of  discs,  which  consists  of  the  unit  disc  Bj,  and  of  k  addi- 
tional small  discs  Bn,  ■  ■  ■  ,Bi;^i  all  of  radius  p,  such  the  centers  of  these  discs  are  placed  on  the 
boundary  of  Bj  at  equally  spaced  positions.  U  p  ia  chosen  to  be  sufficiently  small  (e.g.  of  the 
order  0(~^))  then  it  is  easily  checked  that  for  each  j~2,...,k-i- 1  the  discs  B^  and  B^  are  adja- 

cent  to  each  other  along  the  boundary  ot  the  convex  huli  of  the  iy,  s.  Moreover,  it  can  also  be 
shown  that  for  p  sufficiently  small  each  unbounded  edge  common  to  V'(l)  and  some  other  V(j)  is 
a  full  branch  of  the  corresponding  hyperbola,  and  that  no  two  such  edges  intersect  each  other. 
This  shows  that  Vor{S)  can  have  d,^  iJiiiy  a:  i3(n)  connected  non  r--'^nt3. 

Remark;  The  e.xample  just  given  also  shows  that  the  bouiidary  of  a  single  Voronoi  cell  (V{1)  in 
the  example)  can  have  up  to  0{n)  disjoiiit  connected  cbtn^ohents. 

Suppos:  next  that  for  some  set  S  of  discs  Vor{S)  contain?  a  bounded  component  K.  Then 
the  portion  E  of  a  sufficiently  small  neighborhood  of  K  which  lies  exterior  to  K  must  be  con- 
tained in  a  single  Voronoi  cell  V(i),  since  otherwise  some  arc  of  Vor{S)  would  have  to  enter  any 
such  exterior  neighborhood  of  K,  contradicting  the  assumption  that  A^  is  a  connected  component 
of  Vor{S).  But  if  a  whole  neighborhood  of  /C  in  the  exterior  of  K  lies  in  V(i),  there  must  exist  a 
point  y  6  int(V{i))  such  that  the  line  connecting  y  to  z,  intersects  K,  contradicting  (6).  Thus 
Vor(S)  cannot  have  any  bounded  connected  component. 

Corollary:  Vor(S)  does  not  contain  any  isolated  point.  Moreover,  by  modifying  the  argument 
given  above  one  can  abo  show  that  each  Voronoi  vertex  must  belong  to  at  least  two  edges. 
(Assertion  (11)  below  will  strengthen  this  claim,  by  showing  that  each  such  vertex  must  belong  to 
three  distinct  Voronoi  edges.) 

(10)  No  two  edges  of  Vor{S)  can  be  tangent  to  each  other.  Also,  for  each  point  z  on  a  Voronoi 
edge  e  separating  two  cells  V(i]  and  V(j),  the  segment  connecting  z  to  i,  (or  to  Xj)  is  not 
tangent  to  e. 

Proof:  Since  Voronoi  edges  are  either  straight  segments  or  hyperbolic  arcs  it  follows  that  if  two 
Voronoi  edges  are  tangent  to  each  other  then  any  Voronoi  edges  lying  between  them  at  their 
point  of  tangency  must  be  tangent  to  both  of  them.  Hence  if  there  exist  any  two  tangent  Voronoi 
edges,  then  ihere  exist  two  such  edges  which  belong  to  the  boundary  of  the  same  Voronoi  cell. 
Assume  this  to  be  possible,  and  let  V''(»)  be  a  Voronoi  cell  whose  boundary  contains  two  edges  e, 
e'  which  also  belong  to  V'(;),  V(k)  respectively,  and  which  are  tangent  to  each  other  at  some 
point  y.  It  is  plainly  impossible  for  both  e  and  e'  to  be  straight  arcs;  hence  at  least  one  must  be 
part  of  a  hyperbola.  Let  /  be  the  line  which  is  tangent  at  y  to  both  hyperbolas  (or  to  the  hyper- 
bola and  straight  line)  containing  e,  e'  respectively.  As  is  well  known,  a  line  tangent  to  a  hyper- 
bola at  a  point  y  bisects  the  angle  /ij'/o,  where  /i  and  fn  are  the  two  foci  of  the  hyperbola; 
furthermore,  this  result  abo  holds  trivially  in  case  the  hyperbola  degenerates  into  the  perpendicu- 
lar bbector  of  the  segment  connecting  the  two  points  /i,  f^.  Thus,  in  any  case  I  bbects  the  angle 
between  the  two  segments  connecting  y  to  x,  and  to  x,  (resp.  to  x,  and  Zt).  It  follows  that  the 
these  two  angles  must  be  equal,  and  consequently  the  three  points  j/,  Z;,  Zt  are  colinear,  with  x, 
and  zt  lying  on  the  same  side  of  y.  Suppose  for  definiteness  that  x,  lies  between  y  and  z^.    By 
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definition  of  e  and  e'  we  then  have 

<i(x,,y)-<i(xi,y)=r,-rt 
or 

d{xt,yhd(xj  ,u)=d(xj  ,xt)=rt-r, , 

which  is  to  say,  Bj  is  wholly  contained  in  Bt,  so  that  by  definition  Vor{S)  does  not  contain  any 
edge  bounding  V(j).   This  contradiction  establishes  our  assertion. 

The  second  assertion  follows  from  the  fact  that  no  tangent  to  a  hyperbola  can  pass  through  any  of 
its  foci.  (The  only  exception  is  when  the  hyperbolic  arc  containing  e  degenerates  into  a  half-line, 
but  then  either  B,  or  B,  must  be  whrilly  contained  in  some  other  disc,  so  that  by  convention  e 
does  not  appear  in  our  Voronoi  diagram.)' 

(11)  Let  e  and  e'  be  two  adjacent  edges  along  the  boundary  of  a  Voronoi  cell  V(i).  Then  the 
interior  angle  in  V(i)  between  e  and  e'  is  |«ss  than  180  degrees.  (Stated  otherwise,  in  traversing 
the  boundary  of  V{i)  with  V(i)  to  the  rigt^t,  >ve  make  a  right  turn  as  we  pass  from  one  of  the 
edges  e,  e'  to  the  other.  ,     j      .     .  ^, 

Proof:  We  have  already  shsswB  that->piaad  ici'^qafluat  t«  tangent  to  each  other.  Let  y  be  their 
point  of  intersection,  and  let  /be  the  lii^e  containing  i,,  and  y.  The  two  edges  e  and  e'  cannot 
both  lie  on  the  same  side  of  /  ip  ,the  viciflity.of  y,  bec;?iuse  then  the  segment  connecting  x,  to  a 
point  on  one  of  them  sufficiently  near  y  would  have  to  intersect  the  other  edge,  contradicting  (6). 
Suppose  then  that  e  lies  on  the  kft  side  of  /  (oriented  from  x,  to  y),  and  that  e'  lies  on  the  right 
side  of  /.  Extend  e  along  the  hyperbolic  branch  containing  it  into  the  right  side  of  /  (note  that  a 
hyperbola  is  never  tangent  to  the  segement  connecting  a  point  lying  on  it  to  one  of  its  foci),  and 
denote  the  extended  portion  of  e  by  e"  .  Then  e'  must  lie  between  x,  and  e",  for  otherwise  the 
segment  connecting  x,  to  any  p6int  on  e'  lying  sufficiently  near  to  y  will  intersect  e"  at  a  point  z, 
and  plainly  no  point  z  6  e"  can  belong  to  the  interior  of  V(i)  (because  z  is  equidistant  from  B, 
and  from  some  other  disc  in  5),  again  contradicting  (6).  It  therefore  follows  that  the  interior  angle 
between  e  and  e'  is  less  that  180  degrees,  as  asserted. 

Corollary:  This  argument  shows  that  each  Voronoi  vertex  must  be  incident  to  three  distinct 
Voronoi  edges,  for  if  it  belonged  to  just  two  edges,  at  least  one  of  the  angles  between  these  edges 
would  be  interior  to  some  Voronoi  cell,  and  would  be  greater  than  180  degrees,  contrary  to  what 
we  have  just  shown. 

4.   EfRcient  Construction  of  Generalized  Voronoi  Diagrams. 

Next  we  present  an  algorithm  which  generalizes  the  divide-and-conquer  methods  used  by 
Shamos  [Sh]  and  by  Kirkpatrick  [Ki]  to  construct  the  Voronoi  diagram  of  a  set  of  points,  and 
which  computes  the  Voronoi  diagram  of  a  set  S  of  n  circular  bodies  in  time  0(n  log'n).  The 
algorithm  produces  a  list  of  all  Voronoi  celb  having  nonempty  interior,  and  for  each  such  cell  con- 
structs a  circular  list  containing  the  edges  on  its  boundary,  arranged  in  clockwise  order  (for 
unbounded  celb  this  list  will  also  include  'virtual  edges'  at  infinity  connecting  pairs  of  unbounded 
edges  e,  e'  such  that  the  intersection  of  e'  with  an  arbitrarily  large  circle  lies  immediately  clock- 
wise to  the  intersection  of  this  circle  with  e).  Finally,  the  algorithm  produces  a  table  in  which 
each  edge  points  to  the  two  celb  containing  it. 

For  simplicity,  we  will  assume  in  what  follows  that  no  two  circles  in  S  have  dbtinct  left- 
most points  lying  on  the  same  vertical  line.  If  5  does  not  have  this  property,  then  we  can  apply 
an  infinitesimal  rotation  that  will  make  the  abscbsae  of  the  leftmost  points  of  all  circles  in  S  dis- 
tinct from  each  other.  Cf.  also  [SS]  where  a  similar  technique  based  on  infinitesimal  perturbations 
is  used  to  resolve  degenerate  configurations  arising  in  other  geometric  problems.    The  algorithm 


begins  by  dividing  5  into  two  subsets  R  and  L  of  equal  size,  such  that  the  leftmost  point  w,  of 
each  B,  E  L  lies  to  the  left  of  the  leftmost  point  Wj  of  every  disc  B,  ^  R .  This  partitioning  of  5 
can  easily  be  done  in  time  0{n).  Note  that  it  has  the  property  that  no  body  B,  ^  L  is  wholly 
contained  in  another  body  Bj  E.  R ,  although  a  reverse  containment  is  possible. 

Assume  that  the  Voronoi  diagrams  Vor(R)  and  Vor(L]  have  been  computed  recursively. 
The  main  step  of  the  algorithm  is  to  merge  these  two  diagrams  into  a  single  diagram  Vor{S).  For 
this,  one  must  compute  the  set  C  of  points  y  which  are  simultaneously  nearest  to  a  disc  B,  ^  L 
and  to  a  disc  Bj  E  R .  Following  Kirkpatrick  [Ki]  We  call  C  the  contour  separating  R  and  L  . 
We  will  see  that,  once  C  has  been  computed,  Vor(S)  can  be  obtained  by  taking  the  union  of 
Vor{R),  Vot(L),  and  C,  and  then  by  discarding  (portions  oi)  edges  belonging  to  one  of  the  par- 
tial diagrams  Vor(R),  Vor{L),  whose  points  have  become  nearer  to  some  object  belonging  to  the 
other  set  (these  portions  will  be  delimited  by  the  intersections  ol  these  edges  with  C).  Note  that 
during  this  merging  step  some  Voronoi  cells  V[i)  with  B,  iz  R  may  be  wholly  deleted  from 
Vot(S)  if  B,  happens  to  lie  wholly  inside  sor.-.e  disc  Bj  E  L  . 

Since  C  will  be  a  curve  consisting  of  straight  and  hyperbolic  arcs,  the  complement  of  C  will 
consist  of  finitely  many  open  connected  planar  regions.  Each  such  region  M  is  either  a  union  of 
cells  V(i)  (in  the  finai  diagram  Vor(S))  with  Bi  ^  L  (ifl  which  caie  we  will  call  M  an  L-region), 
or  a  union  of  cells  V(i)  with  S,  6  /?  (in  which  case  M  will  be  called  an  R-region). 

By  assertion  (9)  of  Section  3,  t*  i  coutcur  C  can^  consist  of  several  disjoint  connected  com- 
ponents. Moreover,  it  is  also  possible  for  a  connected  connponent  of  C  to  be  bounded,  as  shown 
by  the  example  appearing  in  Fig.  3. 
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Fig.  3.  A  bounded  connected  component  of  C. 

Fig.  3  shows  that  the  contour  C  can  be  moderately  complex;  the  following  lemma  restricts 
this  complexity,  and  begins  to  develop  some  of  the  facts  that  will  be  needed  to  show  that  tracing 
all  the  components  of  C  need  not  be  expensive. 

Lemma  4.1:  (a)  C  consists  of  a  disjoint  union  of  simple  topologically  closed  curves  without  end- 
points  (i.e.  each  such  curve  is  either  closed  or  stretches  to  infinity  in  both  directions), 
(b)  Let  B,  e  L  and  let  u  be  the  horizontal  half-line  whose  rightmost  endpoint  is  x,.  Then  u  does 
not  intersect  C  (and  consequently  lies  wholly  inside  an  L-regicn). 


(c)  There  exists  precisely  one  L-region;  all  other  componeuts  of  the  compIemeDt  of  C  are  R- 
regions. 

(d)  Each  R-region  has  a  connected  boundary,  consbting  of  a  single  (bounded  or  unbounded)  com- 
ponent of  C . 

Proof:  (a)  It  follows  from  its  definition  that  C  is  a  union  of  Voronoi  edges  of  Vor(S).  It  therefore 
suffices  to  show  that  for  each  Voronoi  vertex  t;  of  Vot(S)  lying  on  C ,  there  are  exactly  two  Voro- 
noi edges  emerging  from  v  which  belong  to  C.  Since  we  have  ruled  out  degenerate 
configurations,  we  can  assume  that  t;  belongs  to  exactly  three  Voronoi  cells  V(i),  V(j)  and  V(k). 
Moreover,  since  v  Q.  C ,  one  of  the  discs  B,,  B,,  B^  must  belong  to  L ,  and  another  of  these  discs 
must  belong  to  R.  Assume  first  that  B,,  Bn  E  L  and  that  B,  Qi  R .  Then,  in  the  neighborhood  of 
V,  the  contour  C  consists  of  the  two  edges  separating  V(j)  from  V(i)  and  V(k)  respectively. 
Much  the  same  argument  applies  it  B,,  B^  ^  R  and  B,  ^  L .  This  proves  (a). 

(b)  Suppose  the  contrary,  and  let  i  be  a  contour  point  lying  on  u.  Let  w  ^  u  be  a  point  at  dis- 
tance r,  from  x,  (that  is,  w  is  the  leftmost  point  of  6,).  It  follows  that  there  exist  discs  Bj  E  R 
(which  is  not  wholly  contained  in  any  other  disc),  and  Bt  6  L  such  that 

rf{z,5/)=  d{z,Bt)<  d(z.B,)--=  d(z,w) 

But  then  some  point  on  B,  (and  in  particular  ita'leftmost  point)  must  lie  to  the  left  of  w,  or  coin- 
cide with  w.  However,  both  these  ^Ssibilities  co&^radict  the  way  in  which  L  and  R  have  been 
defined,  a  contradiction  which  pr6ve3('b).     '  ^     -'    '  ' 

(c)  It  suffices  to  show  that  the  centers  i, ,  Xj  of  each  pair  B,,  Bj  of  dbcs  in  L  are  connected  to 
each  other  via  a  path  which  does  not  intersect  the  contour  C.  Let  u,  (resp.  u^)  be  the  horizontal 
half-line  whose  rightmost  endpoint  is  x,  (resp.  x^),  and  let  w,  6  u,  (resp.  w,  ^  u,)  be  the  leftmost 
point  of  B,  (resp.  Bj).  For  each  /,  let  y,  (resp.  j/^)  be  a  point  on  u,  (resp.  u,)  whose  abscissa  is  t. 
We  claim  that  if  t  is  negative  and  has  a  large  enough  absolute  value,  the  segment  e  =  y^y,  does 
not  intersect  the  contour,  which,  together  with  (b),  implies  that  x,  and  Xj  are  connected  to  each 
other  via  the  polygonal  path  x,y,yjXj  which  is  wholly  contained  in  a  single  L-region.  To  see  this, 
suppose  to  the  contrary  that  there  exists  a  contour  point  z  G  e,  nearest  to  some  B^  E  R  and  to 
some  B„  £  L  .  Then 

j_ 
d(z,B,)  =  d(z,BJ  <  d(z,B,)  =  id"iy,,w,)+  d"'(z,y,))'  ^  d{y,,w,) 

and  similarly 

d(z,Bt)<  d(z.B,)^  d{yj,w,) 

as  l-'oo.  It  follows  that  the  leftmost  point  Wt  of  B^  does  not  lie  to  the  right  of  w,  or  of  w^.  Since, 
by  definition  o{  L ,  R ,  w^  cannot  lie  to  the  left  of  w,,  w,,  it  follows  that  the  three  points  w,,  Wj, 
wic  all  have  the  same  abscissa,  again  contradicting  the  way  in  which  L  and  R  are  defined. 

(d)  Suppose  that  M  is  an  R-region  whose  boundary  consists  of  at  least  two  disjoint  connected 
components  A'i,/vo.  By  part  (a)  of  the  present  lemma  Ki  (resp.  K^)  partitions  the  plane  into  two 
disjoint  components,  both  bounded  by  Ki  (resp.  Ko).  Since  /Cj  and  Ko  are  disjoint,  it  follows  that 
they  both  collectively  partition  the  plane  into  three  components,  one  of  which  contains  M, 
whereas  the  other  two  contain  L-regions.  This,  however,  contradicts  (c),  thus  proving  our  asser- 
tion. Q.E.D. 

In  general  outline,  the  remainder  of  our  argument  is  as  follows.  We  first  show  that,  given 
points  Zf^  on  each  of  the  components  K  of  C,  the  whole  of  C  can  be  traced  in  a  number  of  steps 
bounded  by  0{n).  Next  we  show  how  to  find  such  a  set  of  points  z^-  This  is  done  by  noting  that 
(by  definition)  every  R-region  must  contain  at  least  one  center  z,  of  a  disc  B,  £  R  such  that  B,  is 
not  contained  in  any  other  disc  of  R.  Hence  if  we  iterate  over  all  such  points  x,,  and  connect 
each  one  of  them  by  a  straight  arc  e  to  a  point  of  an  L-region,  these  e  must  together  intersect  all 
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the  components  of  C .  We  will  show  that  in  total  time  0[n  log  n)  we  can  find  such  arcs  e,  each 
intersecting  C  in  just  one  point  z,  and  in  this  way  can  find  %  z  on  each  component  K  oi  C .  This 
leads  to  an  0(n  log  n)  merging  step,  and  hence  to  aa  0{n  log'n)  overall  algorithm. 

To  facilitate  the  tracing  of  C  during  the  merge  phase  of  our  algorithm,  we  will  find  it  con- 
venient (following  [Ki])  to  partition  each  Voronoi  cell  V(»)  (of  either  Vor{L)  or  Vot{R))  into  dub- 
cella  by  connecting  x,  to  each  vertex  v  of  V(i)  by  a  straight  segment  (called  a  spoke,  as  in  [Ki]). 
Clearly  each  subcell  is  an  angular  sector  bounded  by  two  spokes  and  one  Vorouoi  edge.  Note  that, 
given  a  directed  straight  line  or  hyperbolic  arc  e,  its  ini.er"?ct:oa  points  with  the  boundary  of  any 
Voronoi  subcell  can  be  found  in  constant  time,  sssuming  that  we  use  an  appropriate  representap 
tion  of  the  corresponding  diagram  Vor(L)  or  Vor(R). 

Suppose  that  we  have  somehow  found  a  point  z  6  C  (b;it  such  that  z  is  not  in  either 
V'<jr(L)  or  Vor(R)),  for  which  the  two  discs  B,  ^  L ,  Bj  ^  R  nearest  to  z  are  known,  and  sup- 
pose further  that  the  two  subcelb  of  V(i)  in  '/<?r(Z, )  and  of  V(j)  in  Vor{R)  to  which  z  belongs 
are  also  known.  Then  we  can  trace  the  component  K  o.  C  containing  z  as  follows.  We  Erst  find 
the  Voronoi  edge  e  (in  Vor(S))  containing  z.  Note  that  e  is  part  of  the  straight  line  or  hyper- 
bolic arc  equidistant  from  fl,  and  B.,  and  k  aa  edge  lying  on  K.  W.»  begin  tracing  K  by  following 
e  from  z  in  some  direction,  and  by  computing  its  interaectiqn  points  with  the  boundaries  of  the 
two  subcells  U(i),  U{j)  of  V{i),  VI j)  respectively,  coDtaicing,  r .  Su^-rose  for  specificity  that  the 
nearest  of  these  points  along  e  is  the  point  z'  at  which- e  ioterse.'^ts  he  boundary  of  U{i).  If  z' 
lies  on  a  Voronoi  edge,  then  the  contorj-  K  crosses  this  fdge  to  ?nother  Voronoi  cell  V{k)  of 
Vor(L  )  after  z'  (by  assertion  (10)  of  the  previous  section,  two  Voronoi  edges  are  never  tangent  to 
one  another).  In  this  case,  K  continues  after  z'  along  the  Voronoi  edge  e'  containing  points 
equidistant  from  Bj  and  Bj.  On  the  other  hiind,  if  r'  lies  on  a  spoke,  K  wiil  continue  after  ;' 
along  the  edge  e,  but  will  cross  into  another  subcell  of  V{i)  (Note  that  the  contour  can  never  be 
tangent  to  a  Voronoi  spoke,  by  the  second  part  of  assertion  (10)). 

Tracing  the  contour  in  this  way,  we  either  come  back  to  ;■,  in  which  case  the  component  K 
is  a  bounded  component  of  the  contour,  or  else  we  reach  an  unbounded  edge  of  the  contour,  in 
which  case  we  have  to  repeat  the  tracing  procedure  just  outlined  by  starting  from  z  in  the  other 
direction  of  the  edge  e  in  order  to  obtain  the  entire  component  K. 

Let  M  be  the  R-region  bounded  by  K.  Each  cell  V(j)  through  which  K  passes  is  cut  by  K 
into  several  portions,  one  of  which  belongs  to  an  R-region  (and  contains  i, )  while  the  others 
belong  to  the  L-region.  Moreover,  all  the  cell  portions  belonging  to  R-regions  belong  to  the  same 
R-region  M.  Thus,  as  AT  is  being  traced,  we  can  also  note  that  all  these  cell  portions  cut  by  K 
belong  to  M.  Observe  that  A/  may  also  contain  additional  internal  cells  of  Vor(R)  which  have 
not  yet  been  encountered.  These  will  be  dealt  with  during  later  steps  of  the  algorithm. 

Next  we  show  that  if  a  point  Zk  is  available  on  each  of  the  components  K  of  C,  the.  total 
cost  of  constructing  C  is  0{n).  The  complexity  of  the  tracing  procedure  just  described  is  plainly 
0{n['^^+  ni'^^),  where  ni'^^  is  the  number  of  Voronoi  edges  in  K,  and  where  iiA'"^'  is  the  number 
of  intersections  of  K  with  Voronoi  spokes  (in  either  \'or(L)  or  Vor{R)).  As  in  [Ki],  we  can  show 
that  the  sum  over  all  K  of  the  quantities  n}^'+  nl'^'  is  0{n).  Namely  we  have 
Lemma  4.2;  Each  Voronoi  spoke  (in  either  Vor{L)  or  Vor(R))  b  intersected  by  the  contour  C 
in  at  most  one  point.  Moreover  this  remark  also  holds  for  any  segment  one  of  whose  endpoints  is 
the  center  i,  of  some  Voronoi  cell  V{i)  (in  either  diagram)  and  v/hich  is  wholly  contained  in  V(i). 
Proof:  Let  e  be  a  Voronoi  spoke  of  a  cell  V'(t)  in  Vor{L).  Then  each  point  z  at  which  e  and  C 
intersect  each  other  must  lie  on  the  boundary  of  the  cell  V{i)  in  Vor{S),  and  since  this  cell  is 
star-shaped  with  respect  to  x,  it  T  Hows  that  at  most  one  such  intersection  point  can  exist. 
Q.E.D. 

Coi-oUary:  The  total  time  required  to  trace  the  contour,  given  a  point  z  (and  the  two  subcelb  in 
Vor(L),  Vor(R)  containing  it)  on  each  of  its  connected  components,  is  0{n). 


Proof:  The  total  number  of  ed^es  on  C  is  0(n),  because  (7  is  a  subset  of  Vot[S).  The  total 
number  of  intersections  of  C  with  Voronoi  spokes  is  also  0(n),  by  Lemma  4.2.  Hence  the  total 
complexity  of  the  tracing  procedure  applied  to  each  component  /C  of  (7  is 

K 

Q.E.D. 

The  problem  that  now  remains  is  that  6f  finding  a  representative  point  Z/f  on  each  com- 
ponent K  of  the  contour  C  (and  also  finding  the  subcells  containing  Z/f).  For  this,  Kirkpatrick 
[Ki]  uses  a  technique  which  traces  edges  of  mioimum  spanning  trees  for  R  and  L.  However,  this 
technique,  which  works  nicely  for  a  set  S  of  points,  b  not  easily  generalizable  to  sets  of  more  gen- 
eral objects  like  the  circles  which  now  concern  us.  We  will  therefore  present  an  alternative 
approach,  which  works  for  sets  of  circular  d^scs,  but  whose  complexity  is  0(n  log  n),  instead  of 
the  linear  complexity  of  Kirkpatri^k's  technique. 

We  iterate  over  all  the  cells  V(i)  of  Vor(R)  (and  their  corresponding  circles  B,),  proceeding 
as  follows.  Let  B,  E  L  be  such  that  i,  €  ^(y)  (in  Vot{L))  so  that  B,  is  the  circle  of  L  which  is 
'closest'  to  J,.  If  B,  contains  all  o^^Bi,■then,  by  definition,  the  final  diagram  Vot(S)  will  not  con- 
tain a  cell  V(i).  In  this  case^  we  siMply  d'd  not  us«  the  pair  i, ,j^  to  find  a  point  on  the  contour, 
but  go  on  to  consider  the  other  discs  B,i  ER^ 

Next  suppose  that  fl^  does  pot'  tfohtkin' all  of  B,.  li'hen  no  other  disc  B,  ^  L  can  have  this 
property,  and  so  it  follows  that  the  disc  Jj,  in  5  for  wKich  d(x,,Xi]-ri  reaches  its  minimum  belongs 
to  R .  Similarly,  (f(i^,z,)-r,  attains  i(9  minimum  when  /  =  ;,  since  no  disc  in  R  can  contain  the 
whole  of  a  disc  in  L,  and  since  the  recUrsive  construction  of  Vor{L)  described  in  the  following 
paragraphs  will  eliminate  discs  that  are  wholly  contained  in  other  discs.  Thus  the  segment 
e  =  2, 1  must  contain  a  point  z  for  which  d(z,Xi)-ri  reaches  its  minimum  simultaneously  for 
some  Bi  ^  L  and  for  some  other  disc  B,'  ^  R ,  which  is  to  say,  a  point  z  on  C.  Moreover,  since  e 
emanates  from  x^  and  is  wholly  contained  in  V{j),  it  follows  by  Lemma  4.2  that  the  contour  C 
cannot  intersect  e  in  more  than  one  point.  Hence  e  intersects  C  in  exactly  one  point.  Note  that 
the  entire  segment  e  is  contained  in  a  single  subcell  of  V[j)  in  Vor(L  ).  All  that  we  have  to  show 
is  that  we  can  either  find  z,  or  assure  ourselves  that  a  point  on  the  same  component  AT  (of  C)  as 
z  has  already  been  found,  in  total  time  0(n  log  n). 

We  can  proceed  to  find  the  unique  intersection  z  of  e  with  C  using  a  technique  quite  similar 
to  the  contour-tracing  procedure  described  above.  That  is,  we  first  find  the  Voronoi  subcell  (of 
V'{t))  in  Vor(R)  containing  points  on  e  near  x,  (since  e  emerges  from  2,,  this  amounts  to  finding 
the  two  Voronoi  spokes  of  V{i)  between  which  e  lies).  We  then  find  the  intersection  of  e  with  the 
Voronoi  edge  bounding  that  subcell,  beyond  which  e  crosses  into  another  subcell  of  Vor{R).  (In 
the  extreme  case  in  which  e  coincides  with  a  Voronoi  spoke  of  V'(!),  e  will  exit  the  two  subcells 
of  V(j)  in  which  it  lies  at  the  Voronoi  vertex  which  is  the  other  endpoint  of  the  spoke;  it  is  then  a 
bit  more  complicated,  but  still  straightforward,  to  determine  the  Voronoi  subcell  into  which  e 
enters  past  this  vertex.)  Continuing  in  this  manner,  we  partition  e  into  subsegments  e^,  .  .  .  ,e„, 
each  of  which  is  contained  in  some  Voronoi  subcell  of  Vor(R)  or  lies  along  a  spoke  of  Vor(R). 
As  all  this  is  done,  we  keep  track  of  all  the  cells  V{k)  of  Vor(R )  that  have  already  been  encoun- 
tered. If  such  a  cell  is  encountered  for  the  second  time,  tracing  of  the  sequence  of  edges 
ei,  .  .  .  ,e„  stops  immediately  (this  rule  is  justified  by  Lemma  4.3  below).  This  guarantees  that 
the  total  cost  of  traversing  subcells  of  Vor{R)  \s  bounded  by  the  total  number  of  such  subcells, 
and  hence  by  0{n). 

For  each  t  =  l,  .  .  .  ,m  let  V(it)  be  the  cell  in  Vor(R)  containing  e,.  As  tracing  proceeds 
through  the  cell  V{it),  we  check  whether  V(i.)  has  been  encountered  before,  and,  if  not,  whether 
there  exists  z  6  e,  such  that  d(z,Bj)  =  d(z,B,^)  (this  can  be  done  in  constant  time).  As  already 
shown,  there  will  exist  a  unique  point  z  on  e  having  this  property,  and  this  z  will  be  the  required 


intersection  point  of  e  with  C  (the  algorithm  will  reach  this  z  only  if  tracing  is  not  abandoned 
earlier,  because  a  previously  encountered  cell  of  V'or(/?)  is  encountered  again).  Let  e^  be  the  sub- 
segment  of  e  containing  z.  Note  that  z  is  found  in  time  0(p+  d\  where  p  is  the  number  of  sub- 
cells  of  K(i)in  Vot[^R). 

To  show  that  tracing  of  e  can  be  abandoned  as  soon  as  any  cell  of  Vot{R\  is  encountered 
for  the  second  time,  we  will  use  the  following 

Lemma  4.3:  Let  s  be  as  in  the  preceding  paragraph.  Then  for  each  i<i,  Vot[S)  contains  a  cell 
whose  center  is  i,  ,  and  the  point  i,  lies  in  the  same  R-region  as  2, . 

Proof:  Let  A/  be  the  R-region  containing  z^  =  i,,  aid  let  \<t<3.  Pick  any  point  z.  6  ft  (but 
in  case  i  =■  s ,  z^  must  lie  between  z,  and  z\  Since  i<»,  the  segment  z, Zj  does  not  intersect  the 
contour,  and  therefore  is  wholly  contained  in  M.  As  always,  let  B,  be  the  disc  in  L  corresponding 
to  the  cell  V(jt)  of  Vot{R),  and  let  z,  be  its  center.  Suppose  for  the  moment  that  we  have 
already  shown  that  a  cell  '/(t,)  (with  center  z, )  appears  in  Vot{S),  i.e.  that  5,  b  not  wholly 
contained  in  any  disc  of  L  .  The  segment  /  ~  2,2^  is  contained  in  the  cell  V(i,)  of  'Vot[R  ),  which 
is  star-shaped  with  respect  to  its  center  z. ,  Ly  propjrty  (4)  of  Sertioa  3.  Moreover,  since  this  seg- 
ment emanates  from  z,  ,  it  can  intersect  the  contour  in  at  most  one  point,  by  Lemma  4.2.  But 
such  an  intersection  is  impossible,  becauic  faoih  the  endpdiatsof  /  lie  in  tn  R-region  (z,^  lies  in  an 
R-region  because  by  assumption  it  beioB^  to  V(»V)  in  V%t(5)}.  Tbeiffore  the  polygonal  path 
z,  zxx^^  does  not  intersect  the  contour.  But  z,  and  J,  are  coonected  to  e:',oh  other  via  this  path, 
and  hence  lie  in  the  same  R-region. 

It  only  remains  to  show  that  a  cell  V(j()  with  center  z,  appears  in  Vor(5).  For  this,  note  that 
the  point  z^  lies  in  an  R-region,  so  that  it  is  nearer  to  some  disc  of  R  than  to  any  disc  of  L . 
From  this  it  b  plain  that  z,_  is  nearer  to  B,  than  to  any  other  disc  of  5.  But  then  ^,  must  be  an 
interior  point  of  the  cell  V(»()  in  Voro(5),  so  that,  by  definition  of  Vot{S],  a  cell  V[ix)  with 
center  z,   appears  in  this  diagram,  as  asserted.  Q.E.D. 

As  we  apply  the  procedure  just  described  to  each  of  the  discs  B,  ^  R ,  one  of  the  following 
three  situations  will  arise:  either 

(a)  We  discard  B,  immediately,  because  the  nearest  disc  B,  in  L  wholly  contains  B,;  or 

(b)  While  tracing  subcelb  of  Vor(R)  crossed  by  the  segment  e  =  x,Xj  (where  B,  e  L  b  the  disc 
for  which  z,  £  V{j)),  we  encounter  a  subceli  of  some  cell  V'(r)  whose  R-region  M  in  Vor{S)  was 
encountered  before.  In  thb  case  we  conclude  from  Lemma  4.3  that  B, ,  as  well  as  every  other  disc 
of  R  whose  cell  in  Vor{R)  has  been  crossed  by  e  so  far,  lies  in  the  R-region  M.  In  this  case  we 
can  stop  tracing  e  and  go  on  to  process  other  discs  of  R ,  since  we  can  be  sure  that  the  component 
of  the  contour  C  intersected  by  e  has  already  been  explored.  The  algorithm  will  also  note  that 
all  cells  of  Vor{R)  crossed  by  e  so  far  belong  to  M,  to  avoid  repeated  processing  of  these  cells 
later  on.  (Note  that  thb  case  will  arbe  only  when  V(i)  b  an  inner  cell  in  A/,  i.e.  a  cell  not  inter- 
sected by  the  contour);  or 

(c)  The  tracing  procedure  continues  till  an  intersection  r  of  e  with  the  contour  is  found.  In  this 
case  only  new  cells  of  Vor{R)  are  being  traced,  and  z  will  lie  on  a  new  component  of  C.  As  in 
(b),  we  take  note  of  the  fact  that  all  celb  crossed  by  e  during  thb  tracing  belong  to  the  new  R- 
region  just  found,  to  avoid  repeated  processing  of  these  celb  later  on. 

These  observations  imply  that  we  can  find  a  representative  point  on  each  component  of  the  con- 
tour in  total  time  0{n),  provided  that,  for  each  B,  e  R .  the  subceli  of  the  cell  V(»)  of  Vor(L) 
containing  z.  is  already  known. 

To  obtain  this  final  item  of  information,  we  can  use  a  simple  plane-sweeping  algorithm, 
similar  to  those  described  by  [Sh],  [BO],  [f^IP].    The  algorithm  sweeps  the  plane  from  left  to  right 
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and  maintains  a  vertical  "front"  T(a)  consisting  of  the  segments  lying  along  the  line  x  =  a  and 
delimited  by  the  points  of  intersection  of  this  line  with  the  edges  and  spokes  of  Vor(L).  The 
structure  of  T{a)  will  change  only  at  points  a  which  are  either  abscissae  of  Voronoi  vertices  of 
Vor{L),  or  abscissae  of  centers  of  discs  in  L ,  or  points  for  which  the  line  i  =  a  is  tangent  to 
some  Voronoi  edge  of  Vor(L).  The  number  of  such  'transition  points'  is  plainly  0(n),  and  the 
total  number  of  segments  of  T(a)  is  also  0{n)  for  any  real  a.  To  start  the  algorithm,  sort  the 
set  A,  consisting  of  all  transition  points  of  Vor{L)  and  all  centers  x,  of  discs  B,  6  /?,  by  their  x- 
coordinates,  and  initialize  the  list  T(a)  33  a  2-3  tree  for  some  large  enough  negative  real  a.  Both 
these  tasks  can  be  done  in  time  0[n  log  n).  Then  scan  A  from  left  to  right.  For  each  a  £  A,  it 
a  is  a  transition  point  of  Vor{L),  update  the  list  T  by  an  appropriate  combination  of  deletions, 
insertions,  and  merge  operations  applied  to  segments  in  T;  this  can  be  done  in  time  0(k^  log  n), 
where  k^  is  the  number  of  segments  whicji  undergo  these  changes.  Note  that  if  a  is  the  abscissa  of 
a  center  z,  of  some  disc  B,  ^  L,  then  k^  is  the  number  of  Voronoi  edges  on  the  boundary  of  the 
cell  V(i)  in  Vor{L  ),  which  may  be  large.  Nevertheless,  the  total  sum  of  all  the  k^'s  over  all  tran- 
sition points  a  is  always  0(n).  If  a  is  the  abscissa  of  a  center  c  of  some  disc  in  R ,  search  T  to 
find  the  segment  in  T  containing  z,,  from  which  the  Voronoi  subcell  of  Vor(L)  containing  c  is 
readily  obtained.  Proceeding  in  this  way,  we  locate  all  the  centers  of  discs  of  ^  in  Vor{L)  in 
time  0(n  iog  n).  -        ,„    :. 

Together,  the  detaib  just  describi^d  yield  the  following  algorithm  for  constructing  Vor{S): 

1.  Split  5  into  two  equal-size  subsets  L  ,R  aiich  that~the  leftmost  point  of  each  B,  ^  L  lies  to 
the  left  of  the  leftmost  point  of  every  Bj  .6  R  (we  have  assumed  that  no  two  leftmost  points 
have  the  same  abscbsa). 

2.  Compute  V'or(L  )  recursively. 

3.  Apply  the  plane-sweeping  procedure  described  above  to  locate  the  subcell  V{j)  of  Vor(L) 
containing  z,  for  all  centers  z,  of  discs  B,  E  R  ■  Discard  the  disc  B,  £  R  if  it  is  wholly  con- 
tained in  S_,. 

4.  Let  /?'    be  the  remaining  set  of  discs  of  R .    Compute  Vor(R'  )  recursively. 

5.  Construct  the  'contour'  C  as  follows: 

For  each  disc  B,  E  R  whose  R-region  (in  Vor[S))  has  not  yet  been  identified 

a.  Connect  x,  to  the  center  x^  of  the  disc  Bj  S  L  Whose  Voronoi  cell  V(j)  in 
Vor'yL  )  contains  x,.  If  Bj  wholly  contains  B,,  then  Vor(S)  will  not  contain  a  cell 
corresponding  to  B,,  and  we  go  on  to  process  other  discs  of  R . 

b.  Find  the  unique  intersection  :  of  the  contour  with  the  segment  e  =  x.x,  by 
applying  the  tracing  procedure  described  above.  If  that  procedure  detects  an 
intersection  of  e  with  a  cell  V(it)  of  Vor{R)  whose  R-region  A/  has  already  been 
found,  it  assigns  M  as  the  R-region  of  B,  and  of  all  other  discs  of  R  whose  cells 
in  Vor{R)  have  been  crossed  by  e  before  V(k)  has  been  reached,  and  continues 
with  the  main  loop  of  this  phase. 

c.  Trace  the  whole  contour  component  K  containing  2.  An  R-region  indication 
is  thereby  assigned  to  all  discs  B^  £  R  whose  cells  in  Vor(R)  are  encountered. 

6.  Obtain  the  final  diagram  Vor(S)  by  taking  the  union  of  Vor{R)  and  of  Vor(L )  with  C,  and 
then  by  discarding  (portions  of)  edges  of  Vor(R)  or  of  Vor(L)  which  are  cut  off  from  their 
cells  by  C. 

The  algorithm  just  sketched  runs  in  time  0(u  log"n).  Its  costliest  phase  is  step  3,  which 
locates  subcells  of  Vor(L  )  containing  the  centers  of  discs  of  R  . 

Note  that  it  is  a  simple  matter  to  modify  the  algorithm  so  that  it  also  produces  a  mapping 
contain,  which,  for  each  dbc  B,  deleted  from  the  diagram  by  the  algorithm,  gives  the  disc  Bj 
containing  B,  as  found  in  step  5.a  of  the  algorithm. 
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5.   Applications  of  the  Generalized  Voronoi  Diagram;  Possible  Extensions. 

In  this  section  we  show  how  the  generalized  Voronoi  diagram  can  be  used  to  solve  some  of 
the  intersection  problems  mentioned  in  the  introduction  to  this  paper. 

Suppose  that  the  generalized  diagram  Vor{S)  for  a  set  5  of  circular  discs  has  been  con- 
structed, and  that  it  is  represented  by  the  data  structures  described  in  Section  4. 

First  consider  the  problem  of  detecting  the  existence  of  an  intersection  between  any  pair  of 
discs  in  5.   This  can  be  tested  using  the  following  orncedure: 

(a)  First  check  whether  there  exists  a  Voronoi  c?'l  V(i)  having  empty  interior  (i.e.  a  cell  which  is 
missing  from  Vor{S)).  If  so,  B,  is  wholly  cr-ni.aiaed  in  some  other  disc,  and  an  intersection  has 
been  found.  Otherwise,  for  each  edge  e  h-  V:r'S)  belorg'ng  to  the  common  boundary  of  two 
Voronoi  cells  V(j)  and  V(j),  compute  the  value 

p(e)=min{d(x„y\-r,  :  y  €  e}  =  mia{d(Xj,y)-rj  :  y  6  t}. 

If  P(^)^0  fo"^  some  e  6  Vor(S),  then  it  is  clear  that  the  discs  in  5  intersect.  On  the  other  hand, 
if  p{e)>0  for  each  e  6  Vor{S),  then  no  two  discs  in  5  intersect.  Indeed,  suppose  that  two  discs 
B,  and  Bj  intersect  each  other.  Let  /  be  the  segment  (i^j^^^jiForeach  y  €  /  consider  the  function 

f(y)=mh  {d%,y)~xt:  t=5il, .....  .  .n}. 

Note  that  /(2/)<0  for  each  y  E  I,  because  e^cb  such  y 'lies  either  in  B,  or  in  B,,  so  that  either 
d(x,,y)-r,  or  d{Xj,y)-rj  is  <0.  Moreover,  by  ^sumptibD  6oth  V(t)  and  V{j)  have  nonempty 
interiors,  which  implies  by  (2)  that  x,  oelongs  to  V'(i)  arfd  to  no  other  cell.  Hence  /  must  intersect 
Vor(S)  at  least  once.  Let  j/  6  /  be  a  point  belonging  to  some  edge  e  of  Vor{S),  and  let  V(k), 
V(l)  be  the  two  Voronoi  cells  containing  .:  ir  their  boundary.  Then  we  have 

p(e)<d(xt,y)-rt=f(y)<0. 
from  which  our  claim  follows  immediately. 

It  is  easy  to  see  that  a  simple  modification  of  the  procedure  just  outlined  yields  a  solution  to 
the  more  complicated  problem  in  which  the  discs  in  S  are  of  several  colors  and  we  want  to  detect 
intersection  between  two  discs  of  different  colors.  The  appropriate  procedure  in  this  case  is 

(a)  First  check  whether  there  exists  a  disc  5,  which  is  wholly  contained  in  another  disc  B,  of  a 
different  color,  i.e.  if  V{i)  has  empty  interior,  and  x,  belongs  to  a  cell  V{j)  of  a  disc  with  a 
different  color.  Once  the  Voronoi  diagram  has  been  constructed  by  the  method  described  in  the 
preceding  section,  and  has  been  supplemented  by  the  mapping  contain,  we  can  detect  such  cases 
in  0(n)  time.  Note  that  not  every  containment  of  a  dbc  B,  of,  say,  red  color  in  another  disc  B, 
of  a  different  color  can  be  detected  from  the  contain  map,  because  Bj  might  be  contained  in 
another  red  disc  B^,  and  contain  may  map  B,  directly  to  B,.  Nevertheless,  if  B,  is  the  Idrgesl 
possible  disc  contained  in  a  disc  of  s  different  color,  then  contain  will  map  B,  to  some 
differently -colored  disc  containing  it,  so  that  if  any  disc  is  wholly  contained  in  a  disc  of  a 
different  color,  the  procedure  just  described  will  detect  at  least  one  such  situation. 

(b)  If  step  (a)  detects  no  intersection,  compute  the  quantities  p{e),  as  defined  above,  for  all  edges 
e  6  Vor(S).  Then  two  discs  of  different  colors  intersect  each  other  if  and  only  if  there  exists  an 
edge  e  common  to  two  cells  V{i)  and  V{j),  for  which  B.  and  B,  have  distinct  colors,  such  that 
p{e)<0. 

Proof:  Plainly  if  p(e)<0  for  such  an  edge  e,  then  e  must  contain  a  point  y  which  lies  inside 
both  B,  »nd  B,,  so  that  these  two  differently-colored  dbcs  intersect  each  other.  Conversely,  sup- 
pose that  two  differently-colored  discs  B,  and  Bj  intersect  each  other.  Define  the  segment  /  and 
the  function  f  on  it  as  in  the  preceding  paragraphs.  We  can  assume  without  loss  of  generality 
that  V(i)  and  V(j)  have  nonempty  interiors,  for  if  B,  (or  B,)  had  empty  interior  then  it  would 
have  been  wholly  contained  in  another  disc  Bj  of  the  same  color,  so  that  we  could  replace  Bi  (or 
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Bj)  by  Bt  in  what  follows.  Call  the  colors  of  B,,  B,  'red'  ami  'green',  and  call  a  Voronoi  cell  V{p) 
a  'red'  (resp.  'green')  cell  if  Bp  is  colored  red  (resp.  green).  Then  z,  lies  in  (the  interior  of)  a  red 
cell,  whereas  x,  does  not.  Hence  /  must  intersect  Vor{S)  at  an  edge  e  which  separates  a  red  cell 
V{k)  and  a  cell  V(l)  of  a  different  color.  Arguing  as  before,  it  follows  that  p(e)<0  for  this  edge 
Q.E.D. 

Next  consider  the  problem  of  determining  the  shortest  distance  between  any  two  discs  in  5. 
Suppose  that  the  discs  in  S  do  not  intersect  each  other  (if  they  do,  the  above  procedures  will 
detect  this  fact,  and  the  distance  that  we  seek  will  be  0),  and  let  B,,  B,  be  the  two  discs  closest  to 
each  other  among  all  pairs  of  discs  Id  S.  Let  y  fce  the  point  on  the  segment  I=[x,,x  \  equidistant 
from  B,  and  Bj.  We  claim  that  y  6  Vor{S).  For  otherwise,  there  would  exist  another  disc 
Bt  £  S  such  that  d(xt,y)-rt<d{x,,y)-r,:  But  then,  by  the  triangle  inequality, 

d(B,,Bt)=d{x„xt)-rj-rt<d(x,,y)-r,+  d(xt,y)-rt 

<d(x^,y)-rj+d(x„yyr,  =  d{x,,x,}-r,-rj  =  d(B,,Bj), 

contrary  to  assumption.  Thus  y  6  Vor(S).  Moreover,  the  function 

f(z)=mia{d(xt,z}-rt  :  k  =  l,  .  .  .  ,r.} 

attains  its  minimum  value  on  the  wholf  Voronoi  diagram  V'or(5)  at  the  point  y.  This  follows 
since  by  the  triangle  inequality  we  have  2f^{:?)>d{Bii,Bi)  for  each  z  6  Vor(S),  where  V(k)  and 
V(l)  are  the  two  Voronoi  cells  containing  ?,  It  follows  by  the  definition  of  /  and  by  the  preceding 
definition  of  p{e)  that  f{y)  is  the  smallest  of  the  values  p(e),  for  edges  e  of  Vor{S).  Taken 
together,  these  arguments  show  that  t,he  shortest  distance  between  two  discs  in  5  is  equal  to 

min{2^(c)  :  e  an  edge  of  Vor(S)}, 
and  hence  this  distance  can  be  found  in  time  0(n  log^n). 

A  similar  technique  can  be  used  to  find  the  nearest  neighbor  in  S  of  each  B,  G  5.  Indeed,  an 
easy  generalization  of  the  preceding  argument  implies  that  if  Bj  is  the  nearest  neighbor  of  B,, 
then  V(j)  and  V{j)  meet  at  a  common  Voronoi  edge,  and  the  shortest  distance  between  B,  and 
any  other  disc  in  S  is  equal  to 

min{2p(e)  :  e  a  boundary  edge  of  V'(j)}, 
from  which  the  nearest  neighbor  of  B,  is  easily  found. 

These  arguments  extend  easily  to  the  case  in  which  each  of  the  discs  of  5  b  assigned  a  cer- 
tain color,  and  we  want  to  find  the  shortest  distance  between  any  two  diflerenUy-colored  discs. 
For  this,  let  B,  and  B^  be  two  discs  in  5  of  different  colors  such  that  their  distance  is  the  smallest 
among  all  distances  between  two  differently-colored  discs  in  5.  Let  the  colors  of  B,,Bj  be  'red' 
and  'green'  respectively.  Let  y  be  the  point  on  [x,,Xj\  equidistant  from  B,  and  B,.  We  claim  that 
y  6  Vor{S),  for  otherwbe  there  would  exist  another  disc  Bi  such  that 
d{xic,y\-rii<d(x,,y}-r,  =  d{Xj,y)-rj.  If  Bt  is  colored  red,  then  arguing  as  above  we  would  obtain 
d{Bj,Bt)<d{B,,B,).  i.e.  a  shorter  distance  between  a  red  and  a  green  disc.  Similarly,  if  Bt  is 
colored  green,  then  we  would  have  d(B,,Bt)<d{B,,Bj),  again  a  contradiction.  Finally,  if  Bt  b  of 
another  color,  then  both  d{B,,Bt)  and  d(Bj,Bi,)  are  smaller  than  d[B,,Bj),  again  a  contradiction. 

Thus  y  6  Vor(S),  and  similar  arguments  to  those  used  above  imply  that  f(y)  is  the 
minimum  of  all  p(e),  for  edges  e  of  Vor(S)  separating  cells  of  different  colors.  This  shows  that  the 
shortest  dbtance  between  two  differently-colored  dbcs  in  S  is  simply  the  smallest  of  the  values 
2p{e),  taken  over  all  edges  e  of  Vor(S)  separating  two  differently-colored  cells. 

Other  properties  of  generalized  Voronoi  diagrams  deserve  study.  In  particular,  we  would  like 
to  use  these  methods  for  solving  the  second  part  of  Problem  III  mentioned  in  the  introduction 
(that  is,  to  find  the  differently-colored  dbc  nearest  to  any  given  disc  in  S),  and  for  solving  Prob- 
lem   IV,    which    calls    for    an    appropriate    generalization    of   known    point-location    algorithms. 
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